What relaxes muon spins in molecular nanomagnets? 
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We address the cause of the unusual muon spin relaxation (/i + SR) results on molecular nanomag- 
nets (MNMs). Through measurements on protonated and deuterated samples of the MNMs CryMn 
(S = 1) and Crg (S = 0), we show that the muon spin for S 7^ MNMs is relaxed via dynamic 
fluctuations of the electronic spins. A freezing out of dynamic processes occurs on cooling and at 
low temperatures the muon spins are relaxed by the electronic spins which themselves are dephased 
by incoherent nuclear field fluctuations. We observe a transition to a state of static magnetic order 
of the MNM electronic spins in CrrMn below 2 K. 

PACS numbers: 75.50.Xx, 76.75.+i 
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Molecular nanomag nets (MNMs) 0| comprise clusters 
of exchange-coupled transition metal ions which often 
have a negative anisotropy constant favoring a ground 
state with a large eigenvalue of the electronic spin com- 
ponent S z . These systems have been widely studied in 
recent years, most recently in anticipation of their deploy- 
ment as elements of quantum computers 0,3, although 
much interest also centers on the quantum tunnelling of 
the magnetization (QTM) which takes place when the 
magnetic energy levels are at resonance |lfl. The the- 
ory of QTM is based on the Landau-Zener paradigm |4] 
of energy levels brought into resonance by a time de- 
pendent field. The precise details of the mechanism for 
QTM in MNMs remain unresolved and while pairwise 
dipole interactions between spins are certainly important 
[5| a promising suggestion is that resonance is achieved 
as a result of the stochastically varying magnetic field 
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FIG. 1: (Color online.) ^t + SR spectra for protonated and 
deuterated (a) C^Ivln and (b) Cr$ materials, measured at 
T — 4.5 K. Inset: structures of the molecules. 



B n (t), which arises from the incoherent fluctuations of 
nuclear moments surrounding the transition metal ions 
@ (a typical MNM contains ~ 10 2 protons due to the 
organic ligands). 

MNMs have been successfully probed using techniques 
such as magnetization Q , heat capacity Q , neutron scat- 
tering [^, [l(| and electron spin resonance In con- 
trast, measurements made using muon spin relaxation 
(/i + SR) have proven difficult to interpret. Although ini- 
tially it was thought that QTM should be measurable by 



implanting muons into MNMs [ll 13|] , the unambigu- 
ous detection of this effect proved elusive [l4| • Instead, 
/i + SR spectra obtained on high spin systems appeared to 
arise from dynamic fluctuations of a local magnetic field 
distribution at the muon sites, which persisted down to 
dilution refrigerator temperatures [l2, 13, 14 1. Muon re- 



sults on MNM systems all showed similar behavior [14 1 
but it was unclear whether the muon was probing the 
intrinsic behavior of the large electronic spin or some 
residual effect. It was recently argued by Keren et al. 
[HI that ^ + SR is sensitive to the dephasing of the MNM 
electronic spins caused by the incoherent fluctuations of 
nuclear moments in which the metal ions are embedded. 
If this is the case then it makes the muon a valuable 
probe of the potential mechanism behind QTM. In order 
to address the question of what the muon probes in MNM 
systems we have made identical /i + SR measurements on 
protonated and deuterated samples of CryMn (S = 1) 
and Cr 8 (S = 0) [H, E3, El] (structure shown in Fig.H}. 
Here we show (i) that the muon is controlled by the large 
electronic spin in an MNM; (ii) deuteration leads to a 
significant increase in the ^t + SR relaxation rate at low 
temperature in CryMn, implying that muons probe the 
dephasing of large electronic spins by the random mag- 
netic fields due to the nuclei and (iii) that upon cooling, 
a magnetic ground state is reached by a freezing out of 
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dynamic processes that leads to magnetic order in Cr 7 Mn 
below 2 K. 

In a ^ + SR experiment spin-polarized positive 

muons are stopped in a target sample. The time evo- 
lution of the muon spin polarization is probed via the 
positron decay asymmetry function A(t) to which it is 
proportional. Our /i + SR measurements were made on 
the ARGUS instrument at the RIKEN-RAL facility, ISIS, 
Rutherford Appleton Laboratory, UK and on the LTF 
and GPS instruments at the Swiss Muon Source (S/iS), 
Paul Scherrer Institut, CH. Powder samples of the mate- 
rials were packed in Ag foil and mounted on a Ag backing 
plate in 3 He and 4 He cryostats. Typical spectra measured 
for Cr 7 Mn and Crg are shown in Fig. [TJ Above T ps 2 K 
the spectra for all materials differ depending on whether 
protonated or deuterated (Fig. [T]) . 

To compare the Cr 7 Mn and Cr 8 systems, measure- 
ments were made over the temperature range 2 < T < 
100 K. In this regime the spectra for S = 1 Cr 7 Mn 
(Fig. UK a)) were found to be described by the relaxation 
function A(t) = A\ exp(— VAi)+Ab g , where A^g accounts 
for any background contribution from muons that stop 
in the sample holder or cryostat tails. This behavior is 
typical of that observed previously in MNM materials 
llBE! and arises because of the complex dynamic 
distribution of local fields within the material sampled 
by the muon ensemble. The monotonic relaxation and 
the fact that the muons could not be decoupled with an 
applied magnetic field up to 0.6 T places the relaxation 
in the fast-fluctuation limit (20| . 

The spectra measured for the S — Crs samples 
are quite different (Fig D^b)). In this case the relax- 
ation rate is far smaller and resembles a Kubo-Toyabe 
(KT) function /kt(A*) with A = j^(B 2 ) where 7 = 
2tt x 135.5 MHz T _1 is the muon gyromagnetic ratio and 
B is the local magnetic field at a muon site [2(J. However, 
the KT function, which describes relaxation due to static 
random magnetic fields at the muon sites, could not ad- 
equately describe all of the data. This is probably to be 
expected in a complex material such as a MNM where 
many inequivalent classes of muon sites occur and lead 
to a distribution of second moments p(A). The resulting 
muon relaxation is obtained by averaging the KT func- 
tion over this distribution. For simplicity we modeled the 
relaxation by taking p(A) to be a uniform distribution 
up to a maximum Ao- This one-parameter model suc- 
cessfully fitted the data for both materials over the entire 
time range. The fitted values of A are shown in Fig. [2a) 
and are quite T-independent and take average values 
A = 0.56(1) MHz for Cr 8 -h and 0.48(1) MHz for Cr 8 - 
d. Our conclusion from these measurements is that the 
muon is sensitive to the disordered nuclear moments in 
Crs . This is confirmed by the application of a small lon- 
gitudinal magnetic field which quenches the relaxation. 
The larger Ao in Cr§-d compared to Crg-h reflects (albeit 
partially) the larger moment of the deuteron. Most im- 
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FIG. 2: (Color online.) (a) Temperature evolution of the 
relaxation rates A for Cr 7 Mn and Ao for Crg. The lines are 
fits to a phenomenological equation (see main text), (b) Ratio 
of the Cr 7 Mn-h and -d relaxation rates. The line is a guide 
to the eye. 



portantly, the dramatic difference between the measured 
spectra and relaxation rates for S = Crs and S = 1 
Cr 7 Mn samples (Fig. [1} strongly suggests that the muon 
response in MNM systems with S 7^ stems from dy- 
namic fluctuations of the electronic spin. In the absence 
of an electron spin in Crs, the muon spin is relaxed by 
quasistatic disordered nuclear moments. 

We now turn to the role of the nuclei in the spin re- 
laxation process in S = 1 MNMs. The temperature 
dependence of the relaxation rate A for the protonated 
(A h ) and deuterated (A d ) Cr 7 Mn samples is shown in 
Fig. [2^a) . These measurements were made with a mag- 
netic field of 1 mT applied parallel to the initial muon 
spin direction which was intended to quench any residual 
relaxation of the muon spin caused directly by fluctuat- 
ing nuclear spins as seen in the Crs samples. On cooling 
below T ~ 50 K, the relaxation rate A increases before 
saturating below ~ 10 K. with the onset of the increase 
and the saturation occurring at similar values of T for 
both materials. This T-dependence is common to nearly 
all MNM systems that have been previously measured 
with /x+SR 12, ID, 14, 15 1 and is discussed in more de- 
tail below. At high T we see that A d > A h . It is likely 
that at these high temperatures the electronic spins are 
fluctuating very fast and are at least partially motionally 
narrowed from the spectra. Upon cooling the increase 
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in A is greater for the deuterated sample, with A d be- 
coming greater than A h below w 15 K. Most significantly 
the saturation of the relaxation at T < f K occurs with 
A d > A h . The T-dependence of the ratio A d /A h [Fig. Hp))] 
which increases on cooling tends to « 1.7 at the lowest 
temperature. 

The low-temperature results may be understood using 
the model of ^t + SR in MNM systems proposed by Keren 
et al. 15]. Here a muon spin I is coupled to a single 
MNM electronic spin S and dynamic fluctuations of the 
local protons cause S to experience a nuclear stochastic 
field B D (t), along with any externally applied field £? a . 
This system may be described by the resulting Hamilto- 
nian H = g s fi B S ■ [B a + B n {t)} + h<y^I ■ [B a + FS], where 
F couples the muon and electronic spin. We ignore the 
direct response of the muon to the stochastic field, which 
we know to be small from the Crs results. The spin 
relaxation resulting from this model may be simply un- 
derstood. The magnetic field distribution experienced by 
a particular spin population is characterized by a corre- 
lation function whose Fourier transform (known as the 
spectral weight) will be proportional to the spin relax- 
ation rate. We assume that the stochastic magnetic field 
B n (t) experienced by the MNM spin due to the nuclei 
will be described by (B n (t)B n (0)) = (B n 2 ) exp(-t/r n ), 
where r n is the correlation time of the nuclear stochastic 
field. A Fourier transform reveals that the nuclear fields 
will relax the electron spins at a rate 15j 1/t c cx (B^)r n . 



To find the relaxation rate of the muon spin due to the 
fluctuations of the electronic spins, we again assume that 
the correlation function takes the form (£? M (t).B M (0)) = 
(-£?„) exp(— i/r e ), where B^ is the effective local field at 
the muon site due to coupling to the large spin S and r e 
is the correlation time for the electronic spins. Taking a 
further Fourier transform gives the /z + SR relaxation rate 



H El as A = 1 2 M)t c /[1 + ( 7AI |S a |r c ) 2 ] which in the 



limit of small applied magnetic fields becomes 
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The primary effect of swapping protons for deuterons 
might be expected to be the factor 3.26 decrease in the 
size of the nuclear moments experienced by the elec- 
tronic spin and this should be expected to reduce the 
second moment of the field distribution due to the elec- 
tronic spin (B 2 ). In fact there is good evidence for this 
from ESR [3j in the similar system C^Ni, where it was 
found that 1/T 2 decreased by the expected factor of 6 
(corresponding to the reduction in nuclear gyromagnetic 
ratios) upon deuteration, demonstrating that the elec- 
tronic spins are directly relaxed by proton fluctuations. 
Eq. U) shows that a decrease in (B 2 ) will lead to an 
increase in the muon relaxation rate A. The larger mag- 
nitude of A for the CryMn-d sample at low temperatures 
is therefore consistent with the electronic spin being de- 
phased by the nuclei. Since the electronic spins, and 



hence (B^), have the same magnitude in both materi- 
als, then we have A d /A h = (Bl =h )r n=h / (Bl =d )r n=d . We 
note that the measured increase at low temperatures is, 
at most, a factor of w 1.7, rather than a larger factor 
that might be predicted if the distribution width scales 
with the moment sizes alone. It is probable that swap- 
ping protons for deuterons changes not only the local field 
distribution but also its correlation time. 

Finally we address the characteristic temperature de- 
pendence of A observed in all S ^ MNM systems. 
It is generally found in MNMs 14] that on decreasing 
the temperature from T ~ 100 K there is an increase 
in A with the relaxation rate levelling off to some value 
A sa t below ~ 10 K. This behavior may be explained by 
involving two competing dynamic relaxation processes, 
one dominant at high temperatures and described by a 
strongly temperature-dependent correlation time t s (T) 
and one, dominant at low temperatures described by a 
weakly T-dependent correlation time r w . In the pres- 
ence of two competing processes, that with the shorter 
correlation time dominates, giving the smaller relaxation 
rate (since A cx 1/r). At high temperature, therefore, 
we have t s (T) <C t w which results in a strongly T de- 
pendent relaxation which we can crudely model phe- 
nomenologically with A = C exp(U/T). At low temper- 
atures when some of the T-dependent relaxation chan- 
nels have been frozen out, t w <C r s (T) and we have 
A ~ A sat . This behavior results in a phenomenologi- 
cal fitting function 1/A(T) = 1/A sat + l/[Cexp([//T)] 
which has been used previously to characterize these 
MNM systems [Hj]. Fitting this formula to our data 
(Fig. Ufa)) yields C h = 0.23(3) MHz, U h = 46(4) K, 
A^ at = 2.00(4) MHz, C d = 0.070(1) MHz, U d = 51(3) K 
and A d at = 3.3(1) MHz. 

There is a resemblance between the muon results in 
MNM systems and those in some inorganic materials, 
such as CasCoMnOe 2l|, where the physics involves 
significant single-ion anisotropy and a complex (often 
glassy) freezing out of dynamic processes. It is proba- 
ble that the dephasing of the electronic spins we probe 
in MNMs is explainable within the same framework. At 
high temperatures magnetoelastic interactions provide 
the main relaxation mechanism for the electronic de- 
grees of freedom. Our data show U h » U d since, if the 
T-dependence is due to spin-phonon coupling through 
the modulation of local crystal fields, this barrier height 
should only depend on the electronic energy level struc- 
ture. As the temperature is decreased, some relaxation 
channels will be frozen out, increasing the correlation 
time of the electronic moments |22:j . The low tempera- 
ture channel that gives rise to the temperature indepen- 
dent relaxation A sat seen in ^t + SR would appear to be the 
relaxation of the electronic spins by the nuclear fluctu- 
ations, allowing us to identify t w with r e , the electronic 
correlation time discussed above. This situation is simi- 
lar to the electronic T 2 which is dominated by phononic 
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magnetically inequivalent muon sites in the system. The 
results demonstrate that the intermolecular exchange J 
is nonzero in the S = 1 system and, using mean field the- 
ory [25| (assuming z = 6 nearest neighbors) we estimate 
J/k B = 3T c /2zS(S + 1) w 0.2 K. 

This work was carried out at S/iS and the RIKEN-RAL 
facility (ISIS) supported by a beamtime allocation from 
the Science and Technology Facilities Council. We are 
grateful to Alex Amato for experimental assistance and 
the EPSRC (UK) for financial support. 
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FIG. 3: (Color online.) Example spectrum measured at 0.5 K 
for Cr7Mn-h at (S/i + S) showing heavily damped oscillations 
with a frequency of 7 M B/(27r) ~ 30 MHz. Inset: spectra from 
ISIS show a discontinuous change on cooling through below 
2 K. 



contributions at high temperatures and by nuclear con- 
tributions at low temperature [23| . This is not the case 
for the electronic 1 jT\ which is phonon dominated down 
to - 2 K. This difference between 1/Ti and 1/T 2 for ESR 
suggests that the nuclei contribute to the secular part of 
the relaxation, that is, the dephasing of the electron spins 
due to a spread in the net magnetic field at the spin sites 
10. 

In order to further probe the freezing out of the dynam- 
ics in CryMn, measurements were made down to 20 mK 
using the LTF instrument at S/iS. At the lowest temper- 
atures heavily damped oscillations are observed at early 
times (Fig. [3]) in both Ci^Mn-h and -d. Oscillations are 
usually caused in muon spectra by quasistatic magnetic 
order, causing a coherent precession of muon spins. The 
observed oscillations show little temperature dependence 
in the range 0.02 < T < 1 K and are identical for -h and 
-d samples. For measurements made in the same tem- 
perature range at ISIS, oscillations are not discernible 
due to the resolution limit set by the ISIS muon pulse- 
width. However, the spectra do show a discontinuous 
change around T = 2 K [inset Fig. [3], with a loss of 
initial asymmetry and a sharp increase in the apparent 
relaxation rate observed below this temperature. There 
is also significant magnetic hysteresis on the application 
and removal of applied magnetic fields below 2 K. From 
these measurements we estimate that a transition to a 
state of increased static magnetic order takes place at 
T = 1-9(1) K. The heavily damped nature of the oscil- 
lations and the Crg results suggest that there are many 
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